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ABSTRACT 

Looking at the region connecting two clusters is a promising way to identify and study the Warm- 
Hot Intergalactic Medium. Observations show that the spectrum of the bridge between A3556 and 
A3558 has a stronger soft X-ray emission than the nearby region. Suzaku observations could not 
discriminate the origin of the extra emission. In this work we analyze a dedicated Chandra observation 
of the same target to identify point sources and characterize the background emission in the bridge. 
We find that the count number of the point sources is much higher than average field population 
(using CDES 4 Ms as a reference). Moreover, the shape of the cumulative distribution resembles that 
of galaxy distribution suggesting that the point sources are galaxies in a filament. The Suzaku extra 
emission is well explained by the high abundance of point sources identified by Chandra. Eurthermore, 
we used optical/IR observations of point sources in the same field to estimate the density of the 
putative filament as p « 150pt,, below Suzaku sensitivity. 

Subject headings: galaxies: abundances — intergalactic medium — large-scale structure of universe 
— X-rays: diffuse background 


1. INTRODUCTION 

High precisio n cosmological surveys {Planck , ACT, 
and WMAP, see lPlanck Collaborati on et al.|2014[|Sievers| 
et al.||2013 and 


Bennett et al. 


2013 for references) mea¬ 
sured the cosmological baryonic traction as Of, « 0.046. 
In the nearby Universe (z < 2), instead, surveys probe 
only ^ 7 0% of the baryons ( Fukugita et al.|19^ Shull et 


al. 2012|). Large s cale hydrodynamical smiulaLIons (|Cen 
fc Ostriker |19991 [Borgani et al.||2004| [Oppenheimer &: 
Dave [2008 Schaye et al.||2010^ redict that a large frac¬ 
tion of the baryons is in the form of a filamentary gas 
at 10® K< T < 10^ K and density (in terms of baryonic 
density) p < lOOOph. This gas is called the Warm-Hot In¬ 
tergalactic Medium (WHIM) and it is traced by galaxies 
and clusters of galaxies embedded in it. 

The WHIM should be visible in the UV and Soft X- 


ray (^ 100 — 2000 eV) bands . UV (Da nforth k, Shull 
2005 2008 Tripp et al. 2008) and Ly-g ( [Richter et al. 


2004 [Danforth et al. 2010) surveys found evidence of 
several absorbers associated to WHIM at T < 10® K. 


so far there are 

just 

a few dNicastro et al. 2005 Fang 

et al. 120021 20071 Buote et al. 2009 Zappacosta et al. 

20101 JNicastro et al. 

2013D and sometimes controversial 

(Kaastra et al. 

2DDB 

iRasmussen et al.| 

2007) detected 

absorbers. Recently \ 

'Villiams et al. (2013 

suggested that 


the WHIM but rather to intervening galaxies. 

The WHIM is one of several sources that contribute 
to the Diffuse X-ray Background (DXB). Models predict 
that the WHIM contribution is ~ 10% of the total DXB 
(Phillips et al.|2001 Galeazzi et al.|2009 1, making it dif- 
hcult to observe in emission. 

Filaments have been observe d with XMM-Newto n be¬ 
tween clusters A222 and A223 ( Werner et al.|2008 


Diet- 


rich et aL||2012 ), in the visihle\Subaru), X-rays (chan- 


^ corresponding author, galeazzi@physics.miami.edu 


dra) towSds the niassive cluster MACS J0717.5-1-374 5 
( Ebeling et al.|[2004 Ma et al^|2009 Jauzac et al.|2012 ) 
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and with co mbined X-ray/Infrared (ROSAT, Planck) ob¬ 
servations (Planck Collaboration et al. 2013|). All these 
filaments, however, have temperatures ot the order of 
1 keV or more, too high for the WHIM. The analysis of 
the angular correlation function in several “e mpty sky” 
fields provi ded evidence of the WHIM emission (Galeazzi 
et al.||2009), but no detection of individual filaments. 

A region between two nearby cl usters offers the bes t 
chance to detect a filament. Besides Werner et al. (2008), 
there are two reported Suzaku observations towards the 
bridges between A2804 a nd A2811 in the Sculptor super¬ 
cluster ( Sato et ^|2010 ) and between A3556 and A3558 


in the Shapley supercluster (Mitsuishi et al. 2012 M12 
from now on). None of thern found evidence ot WHIM 
emission but were able to set upper limits to the intensity 
of the O VII and O viii emission lines (the best tracers of 
WHIM emission) and to the baryonic overdensity in the 
two regions. 

M12 studied a bridge between clusters A3556 and 
A3558, in the Shapley supercluster, and two nearby tar¬ 
gets, at ^ 1° and ~ 4° away from the bridge, in relatively 
empty regions. The 4°-offset target sets the foreground 
emission (due to Solar Wind Charge Exchange, Local 
Bubble, and Galactic Halo). There is no evidence of 
O VII and O viii emission from the filament and an up¬ 
per limit of 1.5 X 10“^ photons s“^ cm“^ arcmin“^ to the 
O vii intensity, corresponding to a gas at T = 3 x 10® K 
and p ~ 38O(Z/O.lZ0)“^/^(L/3 Mpc)“^/^, was set. The 
spectrum of the on-filament target (see figure shows 
extra emission and a significant enhancement at ener¬ 
gies of Ne IX, possibly indicating an active star forming 
region. The authors attempted to fit the extra emis¬ 
sion with a ^ 2 keV plasma, a temperature much higher 
than the WHIM, but without strong constraints on the 
fitting parameters. The 2 keV plasma model can be as¬ 
sociated to several non-WHIM extended sources (like the 
outskirts of A3556 and A3558 or the plasma in the Su¬ 
percluster itself). 

In our work we investigate the possibility that the 
extra emission comes from unresolved point sources in 
the on-filament target (from now on Shapley filament), 
if they are more abundant than in the reference offset 
targe t s and than in reported literature (Bardelli et al. 
1996 1998). Point sources, in fact, contribute to the 
spectrum with a power law with photon index T ~ 1.4 
(|Mushotzky et al. |2000|). We explore this possibility us- 
ing a dedicated ID ks Chandra A CIS-S observation of the 
on-filament target of M12. The great Chandra angular 
resolution (at the arcsecond scale) makes it possible to 
identify most of the point sources that contribute to the 
spectrum of figure 

We describe the Chandra observation and the reduc¬ 
tion of the dataset in section [21 In section |3] we inves¬ 
tigate the properties of the point sources distribution, 
their contribution to the overall spectrum of the Suzaku 
target, and how they provide limits to the gas density. 

2. THE CHANDRA OBSERVATION 

The Suzaku observation of the Shapley filament was 
performed on July 10 2008 for a total of 30.2 ks of good 
time. The nominal position was (£, 6) = (311.44°, 30.72°) 
and the equivalent absorption density is estimated as 
4.15 X 10^® cm“^. For the analysis of the Suzaku ob¬ 
servation we used the reduced dataset of M12, therefore 



Fig. 1.— The non-Xray background removed spectra and best fit 
models of the off-filament 1° (red circles and line), off-filament 4° 
(blue triangles and line), and on-filament (black squares and line ) 
observations with Suzaku^ as reported in [Mitsuishi et al.| l|2012||. 
The on-filament target shows a strong excess emission tJiat poorly 
fits any thermal models. The extra emission could instead be ex¬ 
plained with excess point sources in the field. The green region 
highlights the Ne IX emission line, particularly strong in the on- 
filament target. The model shown is the “3T” model from M12. 
(A color version of this figure is available in the online journal.) 



Fig. 2.— Suzaku XISl image of the Shapley filament ([0.5- 
2.0] keV). We highli ght the Chandra field of view (magenta square), 
regions removed by [Mitsuishi et'ah ]( |2012[ l (cyan ellipses), and by 
Chandra (green disks). The colorbar scale is in counts s“^. (A 
color version of this figure is available in the online journal.) 

we refer the reader to M12 for a thorough description of 
the data reduction. 

The 10 ks Chandra observation was performed on July 
6 2011 towards {£,b) — (311.36°, 30.67°) using the ACIS- 
S setup to maximize the effective area at energies below 
1 keV. The Chandra field of view covers the innermost re¬ 
gion of the Suzaku pointing, as shown in figurej^ Within 
the Chandra field of view there is only one oithe point 
sources detected by Suzaku, leaving a large area to iden¬ 
tify extra point sources. 

We filtered the dataset to remove bad time intervals 
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plagued by flares. We generated the [0.5 — 2] keV map 
for chipset 7 and the corresponding exposure map (with 
a reference spectrum ba sed on the model described in 
McCammon et al.||2002|). We used the Cl AO wavdetect 
tool to identity the point sources in the map. To reduce 
the risk of false positives in regions with low exposure, we 
masked the maps in order to look for point sources only 
where pixels have exposure higher than 2.5 x 10® cm^ s 
(equivalent to 30% of the maximum exposure), thus re¬ 
ducing the solid angle of the field to ^ 0.017 deg^. The 
detection tool identified 15 X-ray point sources in the 
map, with a detection limit of 5.4 x 10“^® erg cm“^ s“^ 
and a significance threshold for the point-source detec¬ 
tion of 10"®. 

The proper computation of the cumulative count 
numbers requires to estimate the completeness func¬ 
tion. By applying the wavdetect tool to synthetic maps 
of random sources at increasing SB, from 10“^® to 
10“^® erg cm“^ s“^, we obtained the fraction of sources 
detected at a given luminosity. Afterwards, we normal¬ 
ized the count numbers to the completeness function and 
to a solid angle of 1 deg^. 

Figure shows that 4 out of the 15 Chandra point 
sources fall in the mask set for Suzaku point sources re¬ 
moval and one source is exactly at the boundary (we con¬ 
sidered itas being part of the Suzaku mask), therefore we 
extracted the stacked spectra of the 10 previously unre¬ 
solved sources. 


3. THE EXCESS POINT SOURCES 
3.1. Density distribution 

CDFS is the field with the deepest exposure in the X- 
rays and therefore has the best characterization of point 
sources density distribution. We chose CDFS 4 Ms (Xue 


et al. 20111 as a reference to test the cumulative count 
number s distribution of poin t sources in the Shapley fil¬ 
ament. Lehmer et al. (20121 performed a very thorough 
analysis ot the LogN — LogS in the CDFS, dividing the 
point sources in three samples (AGN, stars, and galax¬ 
ies). CDFS 4 Ms gives us the option to compare the 
Shapley filament with the deepest field population avail¬ 
able and to compare the distribution of the detected 
point sources directly with individual populations. 

We tested our method to compute the count numbers 
on the CDF S 4 Ms catalog and compared the results 
directly with Lehmer et al. (2012). Overall our simplified 
approach is i n good agreement w ith the LogN — LogS 
computed in Lehmer et al. (2012), at most our method 
slightly underestimates the number density by a constant 
close to 1 (preserving the shape). 

In figure we show the LogN — LogS distributions of 
point sources in the Shapley filament and in CDFS 4 Ms, 
where the squares represent the data points for the Shap¬ 
ley dataset (black), the total population (red), the AGN 
population (green), and the galaxy population (blue) in 
CDFS, respectively. The Shapley filament clearly ap¬ 
pears to have a much higher population density than the 
CDFS 4 Ms and the slope of the LogN — LogS is possi¬ 
bly steeper than the reference populations. Using the 
number of point sources in the Shapley filament cannot 
be explained with a simple statistical fluctuation, even 
when cosmic variation is taken into account (y^ « 17 
with 6 d.o.f. corresponding to a likelihood of « 9 x 10“®). 



Fig. 3. — Cumulative count numbers of s ources in the Shapley 
filament and in CDFS. For CDFS we used |Lehmer et al.| (|2012j| 
distributions for all point sources, for AGN and galaxies, i ne solid 
line represents the best fit of the Shapley filament LogN — LogS. 
The total count numbers in the Shapley filament is significantly 
higher than the reference CDFS and has a steep slope, similar to 
CDFS galaxies and bright tail AGN. 


Moreover, combined with the lack of excess in the “off- 
filament” targets in M12, it is a strong indication that the 
excess of sources come from the region between A3556 
and A3558 at redshift 0.048. 

Since the data are insufficient to characterize the in¬ 
dividual point sources, we attempted to infer some in¬ 
formation about their nature statistically, once again, 
by a direct comparison with the field distribution in 
CDFS 4 Ms. This ti me we looked at the individual popu¬ 
lations identified by Lehmer et al. (2012|), focusing on the 
galaxy and AGN populations and neglecting stars (the 
fraction of stars over generic point sources is very small 
and the chance of having stars in the Shapley field is 
negligible). We used the slope derived from the fit of the 
LogN-LogS distribution to characterize the source pop¬ 
ulation. Both the Shapley distribution and the galaxy 
distribution wer e fit with a single po wer law. The AGN 
LogN — LogS in Lehmer et al. (|2012 ) is well fit by a bro¬ 
ken power law (and it is a good approximation for the 
distribution of all point sources, too) and we compared 
the Shapley population with the faint and bright tails 
separately. The results of the fits are shown in Table 
The comparison shows that the distribution in the Shap¬ 
ley pointing is incompatible with the faint tail of AGNs, 
but is compatible with both the distribution of galaxies 
and the hard tail of AGNs. However, the luminosity of 
point sources in the Shapley filament is at a much higher 
scale than that actually probed by CDFS. This can be 
explained if we assume that the Shapley point sources be¬ 
long to the region between A3556 and A3558, i.e., they 
are at very small redshift {z « 0.048) compared to CDFS 
(z « I). An excess of galaxies in the filament would, in¬ 
deed, move the distribution in the LogN-LogS plot up 
and right (they would appear brighter). We tested this 
assumption by rescaling the fluxes of all CDFS galax¬ 
ies to the equivalent flux at z = 0.048. The “rescaled” 
LogN — LogS is in good agreement with the Shapley 
sample. On the other hand, it is hard to reconcile our 
observation with the assumption that the excess sources 
are AGNs. Having already discusses how the sources 
must be associated with a low redshift region, the com- 
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TABLE 1 

Deviation betwee n the Shapley f i lamen t distribution and 
THE CDFS 4 Ms|Lehmer et al.|||2012|) in terms of x^, 

P- VALUE, AND (7. 


Sample 


P 

(7 

CDFS (Tot_U/,i„t) 

6.90 

0.01136 

2.47 

CDFS (Tot_Ubright) 

0.91 

0.341 

0.95 

CDFS (Galaxies) 

1.29 

0.255 

1.14 

CDFS (AGN^,i„t) 

7.68 

5.58 X 10-3 

2.77 

CDFS (AGNbright) 

0.25 

0.620 

0.49 


parison with the hard tail of the AGN distribution would 
be inappropriate. The sources in the GDFS are at higher 
redshift, and a direct comparison with Shapley would re¬ 
quire the same rescaling we tested for galaxies. The cor¬ 
rect comparison would therefore be with the faint tail of 
the AGN distribution, but that is incompatible with the 
Shapley distribution. In conclusion, our results indicate 
that the excess emission detected in M12 comes from a 
population of galaxies belonging to the region between 
A3556 and A3558 and could be associated with a WHIM 
filament. 

The galaxy number density in the Shapley filament 
is slightly larger than the number density at the virial 
radii of A3556 and A3558 but the galaxy distribution in 
both clusters is elongated and the density is higher in 
the direction that connects the two clusters (see Fig. 1 
of M12). It has to be noted that the galaxy number 
density at the virial radii of both A3556 and A3558 is 
close to 10 gal Mpc“^ in the NED catalog and that the 
number density of X-ray galaxies in the field (as shown 
in figures and is even larger (~ 70) gal Mpc“^. The 
measured galaxy density is mu ch larger than the average 


densi ty at the cluster radius (Blackburne & Kochanek 


20121 where there should be at most a tew gal Mpc—2 


even assuming that we are at the intersection between 
the virial radii of two clusters. The galaxy density in 
the on-filament target, instead, is remarkably similar to 
the density reported fo r the filament feedin g the cluster 
MAGS J0717.5-b3745 (|Ebeling et al.||2004|), suggesting 
that you are probing a hlamentary structure. 


3.2. Spectral properties 

As confirmation of our results, we verified that the 
spectrum associated with the extra point sources is com¬ 
patible with the extra emission in M12. 

In figure we show the cumulative spectrum of the 
Chandra point sources, the Suzaku “On-filament” excess 
emission model folded through the Chandra response, 
and the cumulative spectra extracted from the CDFS 
point sources. The “On-hlament” emission (red line) is 
much stronger than the CDFS contribution; the emis¬ 
sion from field unresolved sources, therefore, does not 
account for all the extra emission, in particular if we 
consider only the brightest sources (equivalent to a 10 ks 
Chandra observation of CDFS, purple line). The cumu¬ 
lative spectrum of Chandra point sources in the Shapley 
filament (black squares), instead, is in good agreement 
with the excess emission model. This suggests that the 
excess emission is likely due to a number of point sources 
much larger than the average field density. If the point 
sources in this target belonged to a field population, we 
would expect a spectrum in agreement with the sample 



Fig. 4. — Spectrum of the point sources resolved with Chandra 
and of unresolved point sources in the Suzaku field, compared with 
the spectrum of CDFS point sources. Datapoints represent the 
cumulative spectrum of the point sources detected in the Chan¬ 
dra observation. The red curve represents the model spectrum of 
Suzaku unresolved point sources, as reported in M12. The green 
curve represents the cumulative spectrum of point sources in the 
CDFS. The dark and light blue curves represent the cumulative 
spectra of the AGN and galaxies in the CDFS, respectively. The 
purple line represents the cumulative spectrum of all CDFS point 
sources with luminosity corresponding to the detection limit of the 
Chandra observation of the Shapley filament. The contribution of 
the brightest sources in CDFS is too small to explain the extra 
emission in the Shapley filament. All spectra are folded through 
the Chandra response and renormalized to a common field area. 


made of CDFS brightest sources only. 

The fact that the measured point sources spectrum is 
much stronger than a field spectrum suggests that we 
are sampling a non-field population, possibly from a fil¬ 
ament, leaving no room for a contribution from an extra 
thermal component above I keV. 

In the analysis of the filament, M12 found evidence of 
excess emission near the position of the Ne ix emission 
line (see figure [^. Fitting the emission with an extra 
thermal component (attributed to emission from the Su¬ 
percluster) or leaving the Ne ix abundance {vapec with 
the free “Ne” parameter) yielded equally good parame¬ 
ters but there was uncertainty on the detection signifi¬ 
cance due to the uncertainties modeling the broadband 
excess emission. Now we can properly characterize the 
background, modeling it with the extra point sources we 
identified with Chandra, and we have been able to esti¬ 
mate the significance of the Ne ix detection at 2.4cr. Even 
if the Supercluster origin can be ruled out for the extra 
emission beyond 1.2 keV (none of the control Off-targets 
shows sign of it in Fig. [h, we cannot exclude that the 
Supercluster is the cause of the weak Ne ix emission. 

3.3. Limits on overdensity 

In their work, M12 used the Suzaku X-ray observation 
to set upper limits to the gas density toward the line 
of sight (overdensity< 400). In order to estimate the gas 
density in the targeted region, we adopted a different ap¬ 
proach and used data from other wavelengths. Because 
of the poor correspondence between X-ray point sources 
and the NED catalog, we selected all the 0.046 < z < 
0.050 NED galaxies nearby (15’) our field of view to esti¬ 
mate the mass of baryons in galaxies, and from that we 
extrapolated the total baryonic mass. From the catalog, 
we obtained a total of 29 galaxies with photometry in 
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the J and K bands and, in same cases, also in the B and 
V bands (plus 6 galaxies with no photometry). 

We estimated M* in galaxies using the mea s ured m ag- 


nitudes and following Longhetti fc Saracco 


the galaxie s in th e K and V bands and [McGaugh fe| 


(|200^ for 


Schombert (2014) for the B and V bands^ The total 
star mass in the field is M* 1-2 x By adopt¬ 

ing the c losed-box atomic gas-only model of [Baldry et] 


,al. (2008) and aft er a 15% correction to account for un¬ 
resolved galaxies ( Baldry et al.||2008 ) and the 6 galaxies 
with no photometry, we estimated that the baryonic mass 
due to galaxies is ^ 1.4 x lO^^M©. 

In order to estimate the volume sampled, we assumed 
that the NED galaxies belong to a filament. We stud¬ 
ied a section of the hlament extractin g a cylinder of 
depth 3.0 Mpc (Hoshino et al. 2010[ M12) and ra¬ 
dius ^ 0.87 hyg^ Mpc (15’), with depth directed along 
the line of sight. For the cylinder depth, we chose to 
use the size of the diameter of A3556 (the smaller of 
the two clusters). The total volume of the cylinder is 
V « 7.1 Mpc^ « 2.1 X 10®® m® (for the remainder of this 
section all computations already include h dependency). 
This sets the density at pgai « 2.8 x 10“^® kg m“®, equiv¬ 
alent to a baryonic overdenstiy of ~ 70. As a comparison, 
the density measured using weak leasing analysis of the 
filament feeding MA CS J0717.5-1-3745 is p = (206±46)pt, 
( Jauzac et al.|[2M2 |. 

'lb verity that the on-filament target indeed has more 
galaxies than regions away from clusters and filaments, 
we applied the same approach to the two control fields of 
M12 (Offset-1 and Offset-4). Offset-1 has only 1 galaxy 
in the field of view, but no photometry. We used a larger 
NED field of view (20’ of radius), therefore including 
three more galaxies, almost equally spaced from the cen¬ 
ter, in order to obtain upper limits to density. This way 
we computed an overdensity limit of 4pb. In the field of 
view of Offset-4, instead, there are no galaxies so that we 


had to extend the NED field of view up to 30’ of radius 
to find two galaxies and set an upper limit of 2ph. 


According to estimates of the baryon budget (Fukugit a 
|Fukugita & Peebles||2004| |Shull et al.||2012j) 
ike tor ~ 7% of the baryons at low redshifi 


et al. 1998 


galaxies ma^ 


{z < 1), the COM for 5%, and ICM for ~ 4%. Cold gas 
is ~ 1.7%, photoabsorbed (Lya) is ~ 28%, WHIM (OVI 
and Lya) for ^ 30%, and missing baryons are ~ 29%. In 
these works, the galaxy contribution corresponds to the 
stellar mass fraction and the neutral gas well within the 
galactic structure. The CGM, instead, accounts for the 
gas that extends for several hundreds of parsecs around 
the galaxies and, within the large uncertainties of its es¬ 
timate, has the same order of magnitude of the stellar 
mass contribution. In our investigation, we compute the 
halo mass as a fraction (~ 15%) of A4*, much less than 
the ratio between CGM and stellar mass, indicating that 
the gas fraction thus estimated corresponds to the gas 
that more actively takes part into star formation. In our 
attempt to estimate the density in the region targeted 
by the Suzaku observation, therefore, we cannot con¬ 
sider only the stellar mass fraction plus the halo mass, we 
also need to include a separate component that accounts 
for the gas in the CGM and/or in the WHIM closer to 
galaxies and clusters (since at present time there are very 
large uncertainties over the boundaries between these to 
phases). Assuming that the mass in the CGM is the 
same as the mass in galaxies and that there can be ex¬ 
tra contribution from the WHIM, and considering that 
the part of the galaxies are indeed located well within the 
virial radii of A3556 and A3558, it is reasonable to set an 
upper limit to the density of the filament as p « 150pb. 


This work has been supported in part by SAO award 
#G01-12179X. The authors would like to thank R. K. 
Smith and W. Liu for the useful discussion and sugges¬ 
tions. 
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